function [deltaLat, deltaLong] = transform(lat,long)
    x = long -105;
    y = lat -35;
    
    deltaLat = -100+ 2*x + 3*y + 0.2*y.^2 + 0.1.*x.*y + 0.2*sqrt(abs(x))+...
    (2/3)*((20*sin(6*x*pi)+20*sin(2*x*pi))+...
    (20*sin(pi*y) + 40*sin(pi/3*y)) +...
    (160*sin(y./12.*pi) + 320*sin(y.*pi./30)));

    deltaLong = 300 + x + 2*y + 0.1*x.^2 + 0.1*x.*y + 0.1.*sqrt(abs(x))+...
    2/3*((20*sin(6*pi*x) + 20*sin(2*pi*x))+ ...
    (20*sin(pi*x) + 40*sin(pi/3*x))+ ...
    (150*sin(pi/12*x) + 300*sin(pi/30*x)));

end